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ABSTRACT 

We investigate the monopole excitations of the soliton in the Nambu-Jona-Lasinio model. By 
studying the solutions to the corresponding Bethe-Salpeter equation in the background of the 
soliton we exclude the existence of real large amplitude fluctuations. This allows us to treat 
the collective coordinate for the monopole excitations, which parametrizes the extension of the 
soliton, in the harmonic approximation. The canonical quantization of this coordinate yields 
a spectrum which agrees reasonably well with the empirical data for the Roper resonance, 
N(1440), and the corresponding one for the Delta, A(1600). We also comment on going 
beyond the harmonic approximation. 
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1. Introduction 

Bosonization of the Nambu-Jona-Lasinio (NJL) model [1] provides an effective meson 
theory which predicts the properties of the pseudoscalar and vector mesons reasonably well 
[2]. This result has motivated its consideration as a microscopic model for the quark flavor 
dynamics of QCD. Over the past few years it has in addition become evident that the soliton 
approach to the bosonized NJL model suitably describes the properties of the low-lying | + 

and | + baryons [3]. In the present article we will address the question in how far this model 
also permits the description of excited baryons. The natural extension is to consider monopole, 
i.e., radial excitations of the ground states. Such analyses provide descriptions of the Roper 
(1440) and A (1600) resonances, which are excited nucleon and A states, respectively. 

Three different methods have been used to investigate the monopole channel within the 
related (although much simpler) Skyrme model [4, 5]: the scaling approach [6], tt-N phase 
shifts analysis [7], and the linear response theory [8]. In the phase shifts analysis [7] the 
Roper resonance cannot be observed because there is a almost complete cancellation between 
the monopole and rotational channels in the geometrical coupling scheme of ref. [10]. In 
the linear response theory a resonance is observed for the breathing mode at approximately 
400MeV [9]. Unfortunately this resonance cannot immediately be identified with the Roper 
resonance because the coupling to the rotational channel was omitted in that calculation. 
These two approaches appear to suffer from the rotational channel not being treated as a 
large amplitude fluctuation, i.e., like a zero mode. On the other hand the Roper resonance 
is clearly identified in the scaling method, which allows for a dynamical coupling between the 
monopole and rotational degrees of freedom. In this approach the excitation energy of this 
mode comes out at the order of 300MeV [6] which is somewhat too small as compared to the 
experimental value of 500 MeV but nevertheless considerably closer to that value than the 
prediction of the non-relativistic quark model [11]. 

In this paper we will therefore study the scaling method in the framework of the NJL 
model as a first approach to describe excited (non-strange) baryons. We should remark that 
treatments like the phase shift analysis, which involve meson excitations at arbitrary frequen- 
cies, would be plagued by the non-confining character of the NJL model. Once the frequency 
exceeds the binding energy of the valence quarks this quark gets scattered into the contin- 
uum. For commonly adopted parameters of the NJL model the Roper resonance lies above 
this threshold. In the scaling method a collective coordinate is introduced which parametrizes 
the extension of the soliton in addition to those which describe the large amplitude motion of 
the rotational zero mode. The spectrum is then obtained by canonical quantization of these 
coordinates. From the solution to the Bethe-Salpeter equation for monopole excitations of 
the soliton we will argue that an harmonic approximation for the scaling variable is indeed 
justified. In this approximation the feature is circumvented that the NJL soliton is non- 
topological, i.e., by shrinking to zero size the winding number zero sector can continuously be 
reached while the baryon number is carried by three non-interacting valence quarks. 

The present paper is organized as follows. In section 2 we recall the main issues of the 
bosonized NJL model and its soliton solution. In section 3 the solutions to the Bethe-Salpeter 
equation for monopole fluctuations off the soliton are discussed. In section 4 the scaling 
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collective Hamiltonian is determined and quantized. As indicated above the excitation energy 
for the scaling mode is obtained within the harmonic approximation. The numerical results 
as well as related discussions are given in section 5. In section 6 we conclude and comment 
on going beyond the harmonic approximation. The explicit expressions for the kernel of the 
Bethe-Salpeter equation as well as the inertia parameter for the scaling mode are given in 
appendices. 

2. The bosonized NJL model 

The Lagrangian for the NJL model with scalar and pseudoscalar degrees of freedom is 
given as the sum of the free Dirac Lagrangian and a chirally invariant four quark interaction 

[1] 

C = qM-m°)q + 2G £ [(q^-qf + (q^qf) (2.1) 

a=0 \ / 

where q denotes the quark spinor and G is a dimensionful effective coupling constant. m° is 
the current quark mass matrix while the matrices A a /2 represent the generators of the flavor 
group U(Nf). In this paper, we will consider the case of two flavors, Nj = 2, and assume the 
isospin limit m° = m° = m°. Using path integral techniques the model (2.1) can be bosonized 
and expressed in terms of composite meson fields [2]. In Euclidean space with the Euclidean 
time t = it treated as a real number, the resulting effective action is given by 

A = A m + A f . (2.2) 

Here A m and Af are the mesonic mass and fermionic loop contributions, respectively, 



A m = j d 4 x tr (M + M -m (M + M^) + m 2 ), 

A f = TrlogiipE) = Trlog(z^ - (P R M + P L M + j). 



(2.3) 



In eq. (2.3) Pl = (1 — 7 5 )/2 and Pr = (1 + 7 5 )/2 are the usual helicity projectors, while 
"Tr" denotes the functional trace including the traces over color, flavor and Dirac indices. 
Furthermore, Af is a complex matrix which contains the scalar and pseudoscalar meson fields, 
M = S-\-iP . In this work we will neglect fluctuations of the scalar meson field and keep it fixed 
at its vacuum expectation value. In order to determine the minimum of the classical energy 
we will furthermore assume the hedgehog ansatz for the chiral field U . Hence the complex 
field M is given by 

M(r) = mUu(r) = mexp (ir ■ r 0(r)) , (2-4) 

where m = (S) is the constituent quark mass. Demanding the pion decay constant /„. = 
93MeV and mass m v = 135MeV yields the parameters of the NJL model as functions of the 
constituent quark mass [3]. This is a consequence of the gap equation, which determines the 
vacuum expectation value of the scalar fields (S) in the baryon number zero sector of the NJL 
model. We may thus consider m as the only free parameter of the model. 
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Substituting the hedgehog ansatz (2.4) into the expression for the mesonic part A m (2.3) 
of the effective action yields a contribution to the classical energy (subtraction of the reference 
case U = 1 is understood) 

E m = m 2 Jl J d 3 r (f - cos 0(r)) . (2.5) 

Here we have made use of the relation G = m°m/m 2 f 2 obtained from the Bethe-Salpeter 
equation of the bosonized NJL model [2]. In the two flavor model only the real part of 

the fermion determinant (2.3) Ar = ^Trlog( J$ e J/)e) differs from zero. As Ar is ultraviolet 
divergent it must be regularized. We will use Schwinger's proper time regularization [12] which 
introduces an 0(4)-invariant cutoff A 

1 ds + 

A R = -- —Trexp{-sp E p E ). (2.6) 

2 Jl/A 2 s 

For the static meson configurations the energy associated with (2.6) splits into valence quark 
and vacuum contributions [13]. Namely, 

E va \ = ??val|eval| (2.7) 

^vac = -J==Z " • 2.8 

2 ^ Ji/A* V47TS 3 

Here the refer to the eigenstates of the Dirac Hamiltonian 

h = ol ■ p + rnf3 (cos 0(r) + i^fsT ■ r sin 0(r)) , (2-9) 

which commutes with the grand spin operator j + r/2. Furthermore, ?y va i = 0, 1 denotes the 
occupation number of the valence quark, which is the state with the lowest eigenenergy (in 
absolute value). This occupation number has to be adjusted to guarantee unit baryon number, 
i.e., 1 = i] val - (1/2) J2n sgn(e Al ). 

The self-consistent chiral angle s . c .(r) is determined by extremizing the total energy 
functional 1 [14] 

E[0] = E val + E vac + E m . (2.10) 

In the NJL model it turns out that E[Q] depends only moderately on the extension of the 
meson configuration 2 . Whether this gives rise to a large amplitude fluctuation will be subject 
of the next section. 

3. Monopole fluctuations 

In this section we will study the question whether or not the insensitivity of the classical 
energy functional (2.10) with respect to scaling variations of the self-consistent soliton causes 

1 For details of the numerical procedure see ref. [15]. 

2 From figure 6.1 one observes e.g. that E changes by less than 5% when the extension is altered by 30%. 
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the existence of a zero-mode type state in the monopole channel. Such a state would give rise 
to large amplitude fluctuations like e.g. isospin rotations. 

As all our computations are subject to grand spin symmetry the corresponding time de- 
pendent meson fluctuation in the monopole channel is parametrized by [7, 16] 

rj(r,t) = t • rrj(r,t). (3.1) 

Obviously this ansatz describes a pseudoscalar P-wave pion. The expression for the chiral 
field, which contains both the soliton and the monopole fluctuation, reads 

U(r,t) = exp{zr • r (O s . c .(r) + i](r,t))} . (3.2) 

The NJL-model action is next expanded up to quadratic order in the fluctuation rj(r,t). No 
linear term appears because the chiral angle s . c .(r) minimizes the static energy functional. 
As shown in ref. [17] the quadratic term A^ introduces local (/Ci(r)) and bilocal (/C 2 (o;; r')) 
kernels in Fourier space 

A^ = J — I J drr 2 J dr'r' 2 ij(r } cu)/C 2 (o;; r, r')fj(r', — u) 

+ J drr 2 1Ci(r)fj(r\ to)fj(r\ — cu)|. (3-3) 

Here fj(r,u) denotes the Fourier transform of ?y(r,t), i.e., Tj(r,t) = J f^(?", io)e~ tujt . The 
explicit expressions for the kernels of the monopole fluctuation in terms of the eigenvalues and 
-functions of the static Dirac Hamiltonian (2.9) are displayed in appendix A. Here we only 
wish to make a few remarks on the boundary problem. The eigenstates of (2.9) are discretized 
by demanding the upper component of the Dirac spinors ^ to vanish at the boundary (r = D) 
of a spherical box 1 [18]. Eventually we consider D — > 00, in practical computations this means 
that D is significantly larger than the extension of the soliton profile 0(r). The boundary 
condition for ^ transfers to the kernels. As can be observed from the explicit expressions 
given in appendix A the bilocal kernel /C 2 (o;;r, r') vanishes whenever either r or r' equals D 
while K\{r = D) ^ 0. Hence the solutions to the Bethe-Salpeter equation 2 

I dr'r' 2 K,2{w, f'-, r')f]{r\ u) + /Ci(r)?y(r, cu)| = (3-4) 

obey the boundary condition ij(D } uj) = 0. A well defined chiral field U also requires 77(0,0;) = 
0. 

In Fig. 3.1 typical solutions to the Bethe-Salpeter equation (3.4) are shown. In case there 
is no soliton present our solutions are (except of a small vicinity of r = D) identical to spherical 
Bessel functions associated with unit orbital angular momentum 3 . Of course, this is just what 
we expect from a free P-wave fluctuation. 

1 An alternative set of boundary conditions is given in ref. [19]. 
2 In the two flavor case IC'jiu; r, r') is invariant under u> —u>. 

3 This deviation can easily be understood because in the region r k, D the completeness relation 
S/j ^(f)^ (r 1 ) ~ 6(r — r') can only approximately be fulfilled in a finite model space. 
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Electric Monopole 




Figure 3.1: The profiles of solutions to the Bethe-Salpeter equation (3.4) for the constituent 
mass m = 450MeV. These solutions are computed using D = 6fm. The normalization is 
chosen arbitrarily. 
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Two important remarks are in order for the solutions to the Bethe-Salpeter equation (3.4) 
in the background of the self-consistent soliton. First we recognize that the corresponding 
solution to (3.4) deviates from the free P-wave only in the region where the soliton profile is 
non-trivial (r < lfm). Secondly, and more importantly, the resulting energies of the eigen- 
modes in the presence of the soliton agree with those of free P-wave solutions at the order 
of 1%. Usually the existence of large amplitude fluctuations causes a strong reduction of 
the eigenenergy, it may even vanish in the case that the fluctuation corresponds to an exact 
symmetry. We therefore conclude that no large amplitude fluctuation exists in the monopole 
channel despite of the insensitivity of the classical energy with respect to scaling variations. 
Moreover, the spectrum resulting from the quantization of the collective scaling variable will 
be dominated by the properties of the corresponding Hamiltonian in the vicinity of the value 
which minimizes the potential. These considerations justify the harmonic approximation for 
this collective mode which will be discussed in the proceeding section. 

4- The scaling collective Hamiltonian 

The collective breathing mode of the soliton is described by the time-dependent coordi- 
nate x(t) which parametrizes the extension of the meson configuration, x (r, t) = s . c .(x(t)r). 
States with good spin and isospin quantum numbers are generated within the cranking pro- 
cedure which requires collective coordinates R(t) GSU(2) for the (iso) rotations. We therefore 
consider the meson configuration 

M(r,t) = mR^(t)U H (x(t)r)R(t) (4.1) 

where Uu is the hedgehog soliton (2.4). This configuration is substituted into the regularized 
NJL model action and an expansion up to second order in the time derivative is performed. 
From this the collective Lagrangian [6] 

L( X , X) = -E(x) + \a( X )x 2 + \^(x)^ (4-2) 

is extracted. E(x) = -^[®x] * s ^ ne ener gy functional defined in section 2, however, evaluated 
using ® s .c.(x(t) r ) a t a fixed time. I.e., we have substituted the eigenvalues of 

h x = ol ■ p + m/3 (cos X + Z75T • r sin X ) (4-3) 

into eqs. (2.7) and (2.8). The mesonic part of the action (2.3) only contributes to E(x)- The 
inertia parameter a(x) f° r the scaling mode may be interpreted as a position dependent mass 
for the collective coordinate \- Its explicit form in terms of the eigenvalues and -functions of 
h x is displayed in appendix B. Furthermore J? is the time derivative of the collective rotations 

R = (z/2)r • J?. The position dependent moment of inertia a 2 {x) is obtained from the 
expression given in the literature [13, 3] by again substituting the eigenvalues and -functions 
of h x (cf. appendix B). 

In order to avoid ordering ambiguities in the process of quantization we perform a variable 
transformation \ = x(0 such that a(x(0)[^X / d^] 2 = 1- This transformation is allowed as long 
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as a(x) > which actually is found to hold in our numerical studies. It is then straightforward 
to obtain the collective Hamiltonian in terms of the coordinate £ and its conjugate momentum 
P = -id/d£ 



where 



Hj = \P 2 + Vj(i) 



Vj(0 = Vj( X (0) with Vj(x) = E(x) + 



2« 2 (x) 



(4.4) 



(4.5) 



We have already inserted the eigenvalue J of the spin operator, which represents the mo- 
mentum conjugate to R. This procedure is justified because the generators for rotations and 
scalings commute. 

The above motivated harmonic expansion is carried out separately in each spin channel. 
For this purpose the minimum of the potential Vj(£) is determined. We denote its position 
by <^™ n . Then the spectrum of the collective Hamiltonian (4.4) is approximated by 



E, 



VAC. 



min\ 
J ) 



+ ujj [n + 



1 



\ 



dWj(0 




vnxf 11 ) 






a( X f n ) 



(4.6) 



with Xj" n = x(Cj un ) an( i the prime indicating a derivative with respect to x- It should 
be remarked that Xj 1 ™ minimizes Vj(x)- The harmonic approximation has obviously the 
advantage that we do not have to explicitly carry out the change of variables x — > £• A 
further justification of this approximation is provided by the fact that our computation of 



a 



x) cannot be generalized to arbitrary large x, cf. the discussion at the end of appendix B. 
The nucleon corresponds to the quantum numbers J = 1/2 and n = while the Roper 
resonance is associated with J = 1/2 and n = 1. As an illustrative example we also list the 
expression for the A-nucleon mass difference 



Ma — M 



1 



N 



(4.7) 



Obviously this mass difference acquires an additional contribution, which is due the treatment 
of the scaling mode quantum variable. 

5. Numerical results 

In this section we present our predictions obtained in the NJL soliton model for the masses 
of excited baryons in the monopole channel using the harmonic approximation to the breath- 
ing mode. It is known that the absolute masses of solitons acquire substantial reductions 
associated with meson loop corrections [20]. Although the absolute mass of the soliton in the 
NJL model is not as large as in the Skyrme model (for parameters fitted to the meson sector) 
these reductions are relevant in the NJL model as well [21]. As they effect all baryon states 
approximately equally we will concentrate on mass differences only. 
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As already mentioned above the constituent quark mass m is the only undetermined pa- 
rameter. Although for a unit baryon number soliton solutions exist for m > 325MeV [14] 
these solutions only represent local minima of the energy E (2.10). Unless m > 420MeV 
the configuration consisting of three non-interacting valence quarks is energetically favored 
against the soliton configuration. As the NJL soliton is of non-topological character it can 
continuously be deformed from the local to the global minimum by shrinking it to zero size 
(x — *■ oo) without encountering an infinite energy barrier. We therefore consider only the 
region in parameter space where the soliton represents the global minimum of the energy 
functional i.e., m > 450MeV. A value of m that large is also mandatory to find a pronounced 
minimum of Vj (4.5) in the A channel (J = 3/2). 

In table 5.1 the numerical results for quantities appearing in the mass formula (4.6) are 
shown for various values of m. We have first computed the position Xj" n °f the minimum of 
Vj(x) (4.5). As the moment of inertia cancels from eq. (4.6) for J = the result Xo 1 ™ = 1 
confirms that s . c .(r) indeed minimizes the classical energy. We observe that the deviation of 
X™ 11 from unity is significantly smaller in the NJL model than the Skyrme model calculations 
[6] which yield X3/2 as smai l as 0.4. Although the breathing mode potential is shallow it is at 
least steeper than in the Skyrme model. Subsequently Xj™ nas been employed to evaluate the 
classical mass E, the moment of inertia a 2 and the breathing frequency uj defined in eq. (4.6). 
Obviously represents our prediction for the mass difference between the Roper resonance 
and the nucleon. From table 5.1 we also deduce that the classical energy E is indeed quite 
insensitive to variations in the scaling variable. Glancing e.g. at the case m=500MeV shows 
that E changes by only about 5% when x is reduced by 20%. On the other hand the moment 
of inertia a 2 and the breathing frequency uj crucially depend on x- 

Now we come to the central issue of this paper, the spectrum of the scaling mode. If 
this mode were not treated quantum mechanically the A nucleon mass difference would be 
equal 3 /2a 2 . The experimental value of this mass difference (293MeV) corresponds to a 2 
5.12GeV~ which is obtained for m ~ 420MeV. However, as already indicated at the end 
of section 4 the treatment of the scaling mode as a quantum variable drastically alters this 
result. From table 5.2 we observe that this mass difference is best reproduced for values of 
the constituent quark mass as large 550MeV. This result also gives an a posteriori 

justification for the harmonic approximation which requires a pronounced minimum of Vj 
(4.5). This is not the case for small constituent quark masses. 

We also find that the parameter m ~ 550MeV not only correctly reproduces the A- 
nucleon splitting but also leads to a reasonable prediction for the mass difference between the 
Roper resonance and the nucleon. As a matter of fact the Roper nucleon splitting is almost 
independent of the constituent quark mass. It appears to be a common feature of the scaling 
approach to soliton models that the Roper comes out on the low side [6]. For m = 550MeV 
the agreement for the first excitation above the A is equally good although its position is more 
dependent on m. In table 5.2 we have also displayed our predictions for second excited states 
in the J = 1/2 and J = 3/2 channels. Also in these cases our results compare reasonably well 
with the experimental data although these are not exactly identified 1 . 

1 In ref. [22] these states are listed as three star resonances. 
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Table 5.1: The parameters of mass formula (4.6) at Xj 1 ™ which minimizes the potential Vj(x) 
(4.5). 



m=450MeV 



J 


^ .nun 
Xj 


£(GeV) a 2 (l/GeV) 


cuj(GeV) 





1.00 


1.232 4.79 


0.450 


1/2 


0.99 


1.233 4.80 


0.400 


3/2 


0.83 


1.269 5.47 


0.266 


m=500MeV 


J 


^ .mm 
Xj 


£(GeV) a 2 (l/GeV) 


cuj(GeV) 





1.00 


1.221 4.17 


0.456 


1/2 


0.97 


1.222 4.24 


0.403 


3/2 


0.79 


1.287 5.46 


0.303 


m=550MeV 


J 


A .mm 
Xj 


£(GeV) a 2 (l/GeV) 


cuj(GeV) 





1.00 


1.208 3.75 


0.461 


1/2 


0.96 


1.210 3.89 


0.405 


3/2 


0.77 


1.294 5.46 


0.324 


m=600MeV 


J 


A .mm 
Xj 


£(GeV) a 2 (l/GeV) 


cuj(GeV) 





1.00 


1.193 3.46 


0.465 


1/2 


0.95 


1.196 3.65 


0.404 


3/2 


0.75 


1.294 5.45 


0.354 



Table 5.2: The predictions of the masses of the low-lying baryons according to the mass 
formula (4.6). Given are the mass differences to the nucleon J = 1/2, n = 0. The energy 
dimension is MeV. 





m 




J 


n 




450 


500 


550 


600 


Expt. [22] 


3/2 





A(1232) 


234 


270 


292 


315 


293 


1/2 


1 


iV(1440) 


400 


403 


405 


404 


501 


3/2 


1 


A(1600) 


500 


573 


616 


669 


661 


1/2 


2 


iV(1710) 


800 


806 


810 


808 


742 - 802 


3/2 


2 


A(1920) 


766 


876 


940 


1023 


962 - 1032 
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6. Conclusions 

We have investigated the monopole excitations of baryons within the NJL chiral soliton 
model. Although the surface of the classical energy is quite flat in the direction of scaling 
variations of the static soliton it has turned out that no large amplitude fluctuation exists for 
the scaling mode. In account of this result we have argued that the spectrum in the monopole 
channel is characterized by the properties which the potential exhibits at its minimum. When 
elevating the scaling mode to a quantum variable we have therefore treated this potential 
in the harmonic approximation. Adopting the constituent quark mass m = 550MeV within 
this approach we have obtained a reasonable agreement with the available data for the mass 
differences of the exited non-strange baryons in the J = 1/2 and J = 3/2 channels. Although 
this value for the constituent quark mass appears to be somewhat high there is nothing special 
about, it just represents the only free parameter of the model. On the contrary, such a 
large value is appreciated since it makes the harmonic approximation more reliable since the 
minimum of the potential is more pronounced thereby providing an a posteriori justification 
of the method. Let us nevertheless comment on treating the Hamiltonian (4.5) beyond the 
harmonic approximation. For this purpose we have displayed the complete dependence of the 
potential Vj(x) on the scaling variable x m figure 6.1. Obviously the potential stays finite 
as x ~ > °°? rather E[Q] — > 3m. This limit just corresponds to the absence of the soliton 
(shrunk to a point) and the baryon number carried by three non-interacting valence quarks. 
Of course, in a confining model such a minimum would not exist. Stated otherwise, our 
harmonic approximation represents a (crude) way to imitate confinement. For our preferred 
value m = 550MeV of the constituent quark mass we find that decays of the first excited states 
(n = 1) into three free quarks are on the border of being energetically forbidden. The higher 
excited states may, after passing through a finite energy barrier, decay. As the full potential 
is more shallow than the one approximated harmonically the predictions for the masses of the 
excited states will be reduced. It should also be mentioned that any treatment which goes 
beyond the harmonic approximation suffers from ordering ambiguities when quantizing the 
breathing coordinate especially because a(x) may contain large derivatives. 

There is an additional feature we can read off from figure 6.1. Although we observe quite a 
pronounced minimum of the potential for the soliton, it becomes the more shallow the larger 
the spin quantum number is. For J > 3/2 a minimum ceases to exist. This fact may be 
considered as an indication that the NJL soliton model does not contain baryons with spin 
larger than J = 3/2 when the rotational degrees of freedom are treated beyond the cranking 
approximation. Of course, this is expected [23] within a model of baryons which is formulated 
in terms of quark degrees of freedom when Nc = 3 is adopted for the number of colors. 

Appendix A: Kernels for monopole fluctuations 

In this appendix we will present the explicit expressions for the kernels which enter the 
quadratic form for the monopole fluctuation (3.3). The expressions quoted in this appendix 
refer to Minkowski space. 

It is suitable to introduce the chiral rotation 

0(r) . A . 0(r) 
T(r) = cos «75T • fsin (A.l) 
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J=0 

J=1/2 

J=3/2 




1.0 



1.5 



Figure 6.1: The potential of the collective Hamiltonian (4.5) for the spin quantum numbers 
J = 0, 1/2 and 3/2 as function of the scaling coordinate \- The constituent masses adopted 
here are m = 450MeV (left panel) and m = 550MeV (right panel). 
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since it permits a compact notation for the Dirac operator in the presence of the monopole 
fluctuations (3.1) 

i/3p= idt-a-p- mT(r)/3 jl + ij 5 t ■ rr]{r,t) - ^(r,t) 2 + . . .Jrt(r). (A. 2) 

When computing the functional trace the chiral rotation can straightforwardly be imposed 
onto the eigenfunctions of the static Dirac Hamiltonian (2.9). This procedure simplifies the 
expressions for the direct coupling between the soliton and the meson fluctuation. Since the 
eigenfunctions of the static Dirac Hamiltonian depend on the soliton profile functionally there 
is also an indirect coupling. 

Substitution of the expansion (A. 2) into the general expressions for given in ref. [17] 
yields the local kernel which also gains a contribution form the mesonic part of the action 
(2.3) 



fCi(r) = — 27rm 7r / 7r cos0(r) — mNcTj 



val 



^vii(r)T(r)/3Tt(r)Vw(r) 



47T 



+ ^miV^ V ai£sgn( e ,)erfc(N) J _^t( r ) T (r)/3Tt(r)^(r), (A.3) 

where we have indicated the average over the angular degrees of freedom. Similarly the bilocal 
kernel is obtained as 

xV> v t al (r)T(r)/3 75 r • rT\r)iUr)4{r')T{r')fj l5 T • r rV)Vw(r') 
+ r *£ %£ ,;,')/ t / ¥ (A.4) 

x,f t(r)T(r),8 T5 r • rT*(r)Vv(r),f 1(t')T{t')0 1s t ■ r'rt(r')^(r'). 

Both, the local as well as the bilocal kernels, are decomposed into valence (~ ?y va i) and vacuum 
contributions according to eqs. (2.7,2.8). As at large |r| the chiral rotation T equals unity 
upper and lower components get connected. This causes £2(^5 r') to vanish when either r 
or r' lies on the boundary. The regulator function in (A.4) is given by a Feynman parameter 
integral 

f°° f~s f e _ s£ ^ — e _ S£lJ 
R(e„e„;u 2 ) = ldsJ — \ (A.5) 

Jl/A 2 V 47T [ S 

+ [cu 2 - (e^ + £j,) 2 ] dx exp [-5 (xe 2 + (1 - x)e 2 - x(l - x)u 2 ^ j 
which describes the quark loop in the background of the static soliton. 
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Although we have left the number of colors Nc as a free parameter, it is always implied 
that it assumes its physical value Nc = 3. 

Appendix B: Inertia parameters 

In order to extract the collective mass a(x) f° r the scaling variable xif) within the NJL 
model we parametrize the time dependence of xif) as 

x(t) = xo + S(t) (B.l) 

and expand the action up to quadratic order in S(t) while keeping all orders in Xo, which 
is assumed to be time independent. There are various ways to perform this computation. 
One might e.g. straightforwardly adopt the treatment of ref. [17] where time dependent 
fluctuations off the chiral soliton have been considered. In the present case a more elegant 
way is to define the translation operator 

f = exp(6(t)d Xo ). (B.2) 

Then the meson configuration (4.1) may be written as 

M(r,t) = mltf(t)f- 1 UH(xor)fR(t). (B.3) 

This parametrization may be transferred to the Dirac operator (2.3) defining the Dirac Hamil- 
tonian 

z/3|)=r- 1 i?t (id t -h x )RT with h x = h Xo + • fl + i6d Xo . (B.4) 

In addition to the Coriolis term an expression involving the time derivative of the scaling 
coordinate has been induced, 8 = X- Since the eigenstates of h Xo are properly normalized the 
translation operator T is actually unitary. This allows us to absorb both operators T and R 
by redefining the quark fields q' = TRq. The collective mass a(xo) is then extracted from 
the term, which is proportional to S 2 = % 2 , of the NJL model action. This computation is 
completely analogous to standard determination of the moment of inertia a 2 (xo) [13] resulting 
in 

a (Xo) = Oval(Xo) + Ovac(Xo) (B.5) 

with 

9 AT sr KHd X o|val)| 2 

2??val/Vc 2^ ~ 

NcJ2\WM\ 2 fe(^e^A) (B.6) 



Oval(Xo) = 
Ovac(Xo) = 
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with the cut-off function 



/©(e^, e„; A) 



2 e -(^M) 2 



sgn(e J/ )erfc 



sgn(e At )erfc 



7T 



e 2 - e 2 



2(e M - e p 



;b.7) 



Actually /©(e^, t v \ A) is proportional to dR{e ll ,e v ]Lo 2 ) / du 2 



This just reflects the fact 



that we have treated x as a small fluctuation and expanded up quadratic order in the time 
derivative. It is important to note that in eq. (B.6) the eigenstates and -values of h Xo have 
to be substituted. In practice it is convenient to employ the identity 



1 



m 



-(H [^xo Ao] \ v ) 

-(fi\/3rQ' s c (xor) [sin s . c .(XoO - i^r • r cos s . c .(XoO] W) 



where the prime indicates a derivative with respect to the argument. The regulator function 
(B.7) vanishes for identical energies. Therefore the limit 



lim 



1 



r/e (e M , e^; A) 



e, 



is finite. 

For completeness we also give the explicit expression for the moment of inertia [13] 



« 2 (Xo) = «val(Xo) + «vac(Xo) 



;b. 9 ) 



Wl 



ith 



«val(X0j 
«vac(Xo) 



1 at KH^MI 2 



(B.10) 



where r 3 denotes a Pauli matrix. Again, \fj) and refer to the eigenstates and eigenvalues of 
h Xo , respectively. One word of caution is necessary when considering the limit Xo ~ > 00 • Then 
the soliton actually is absent and the eigenstates of h Xo are also eigenstates of the isospin oper- 
ator r 3 . Since the diagonal elements of the regulator function vanish [13] /©(e^, e A*! A) = the 
moment of inertia vanishes in the limit Xo ~ > 00 • ur numerical calculation does not exhibit 
this feature. The reason being that states employed to diagonalize h Xo are no eigenstates of 
r 3 rather they are eigenstates of the so-called grand spin operator, which is the sum of the 
total spin and isospin. The perturbation expansion of h x (B.4) cannot straightforwardly be 
generalized to Xo ~ > 00 m the grand spin basis because states with different grand spin eigen- 
value (but identical spin and isopin) become degenerate. Since our techniques are unable to 
remove this degeneracy the small energy denominators in eq. (B.10) cause a 2 to increase for 
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large Xo rather than to decrease. Similar arguments hold for Xo ~ > since h Xo=0 is also isospin 
invariant. These deficiences are fortunately circumvented by the harmonic approximation 
(4.6). 
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